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Abstract: The Tibet air siiower array, which has an effective area of 37,000 square meters and is located 
at 4300 m in altitude, has been observing air showers induced by cosmic rays with energies above a few 
TeV. We have a plan to add a large muon detector array to it for the purpose of increasing its sensitivity to 
cosmic gamma rays in the 100 TeV energy region by discriminating them from cosmic-ray hadrons. We 
have deduced the attainable sensitivity of the muon detector array using our Monte Carlo simulation. We 
report here on the detailed procedure of our Monte Carlo simulation. 
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Introduction 



Supernova remnants (SNRs) are the best candi- 
dates for acceleration of hadronic cosmic rays up 
to the knee in the cosmic-ray energy spectrum at 
^ lO^^eV. Consequently, gamma rays in the 100 
TeV region (10 - 1000 TeV) originating in 7r° de- 
cay following inelastic collisions between acceler- 
ated charged cosmic rays and the ambient medium 
are naturally expected. Since the expected flux of 
7r° decay gamma rays is low, an apparatus with a 
high duty-cycle and a large field of view, e.g. an 
air shower array, is suitable for this purpose. Air 
shower arrays thus far, however, have not been sen- 
sitive enough for the detection of such gamma rays. 
This is because their angular resolution and/or their 
power of discriminating gamma-ray induced air 
showers from hadron-induced ones were insuffi- 
cient. 

The Tibet air shower (AS) array has succeeded in 
detecting celestial TeV gamma rays [1] and is op- 
erating with effective area of 37,000 m^. At the 
100 TeV energy, its angular resolution and energy 
resolution are estimated to be 0.2° and 40%, re- 
spectively. We are now planning to build a large 
water Cherenkov muon detector (MD) array in the 
underground of the AS array for the purpose of 
gamma-hadron discrimination. This future plan 
would be able to open the 100 TeV energy window 
in gannma-ray astronomy. 

Muon Detector Design 

The currently proposed configuration of the MD 
array is shown in Figure 1. It is composed of 12 
pools, each of which consists of 16 cells. Each 
cell is a waterproof concrete tank which is 7.2 m 
wide X 7.2 m long x 1.5 m deep in size. Two 
20 inch-in-diameter photomultiplier tubes (PMTs, 
Hamamatsu R3600) are put on its ceiling, facing 
downwards. Its inside is painted with white epoxy 
resin to waterproof and to efficiently reflect wa- 
ter Cherenkov light, which is then collected with 
the PMTs. The MD array is set up 2.5 m under- 
ground (2.0 m soil + 0.5 m concrete ceiUng, ~19 
radiation lengths) in order to detect the penetrat- 
ing muon component of air showers, suppressing 
the electromagnetic one. Its total effective area 



Figure 1: Schematic view of the Tibet AS+MD ar- 
ray. Open squares and open circles represent the 
surface scintillation detectors that compose the Ti- 
bet AS array. Note that the AS array drawn here 
is upgraded from the current version so that its ef- 
fective area becomes 50,000 m^ by modifying the 
configuration of the scintiUation detectors. Filled 
squares show the proposed Tibet MD array 2.5 m 
underground. 

amounts to 9,950 m^ for muon detection with the 
energy threshold of approximately 1 GeV. The ad- 
vantages of using the water Cherenkov type de- 
tector are high cost performance and its capabil- 
ity to exclude the influence of the environmental 
background radioactivity thanks to the Cherenkov 
threshold. 



Simulation 

The air shower events induced by primary cosmic 
rays and ganrnia rays were generated in the energy 
range from 0.3 TeV to 10 PeV and within zenith 
angle less than 60° along the Crab's orbit using the 
Corsika Ver.6.204 code [3]. We used QGSJETOlc 
for the hadronic interaction model and adopted a 
chemical composition model [4] based on direct 
observational data for the primary cosmic rays. We 
assumed a differential energy spectrum E~'^-^ for 
the primary gamma rays. Air shower events were 
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Figure 2: The light attenuation length (left) and the 
quantum efficiency of the PMTs (right). Both are 
dependent on wavelength of light. 



uniformly thrown within 300 m from the array cen- 
ter This radius is sufficient to collect all the air 
shower events which actually trigger the AS ar- 
ray. The response of the AS array was simulated 
using the Epics uvS.OO code [5], which had been 
already established. From the simulation, we de- 
duced the air shower direction, core position and 
the sum of the number of particles per m^ detected 
in each scintillation counter (^ p) etc. Note that in 
this simulation we used an upgraded version of the 
AS array, the effective area of which was increased 
from 37,000 m^ to 50,000 m^ by modifying the 
configuration of the current scintillation detectors. 
The total number of the scintillation detectors were 
not changed. 

We selected air shower events based on several 
conditions, i.e. software trigger condition of any 
fourfold coincidence in the FT counters recording 
more than 1.25 particle in charge, air shower core 
position located in the array and the residual error 
of direction reconstruction less than 1 .0 m. 

The secondary particles of the surviving air shower 
events were subsequently fed into the simulation 
of the soil absorber, then the MD array. Their 
responses were simulated based on GEANT4 8.0 
code [6]. We assumed that the soil was a mixture 
of 70% SiOz, 20% AI2O3 and 10% CaO with den- 
sity 2.0 g/cm^ and thickness 2.0 m. For the MD 
array, its detailed structure made of concrete (2.3 
g/cm"^, 100% Si02) was taken into account and the 
ceiling thickness of each cell was set to be 0.5 m. 
The reflectance on the inner surface of the cells was 
assumed to be 70% with isotropic reflection. The 
light attenuation length and the quantum efficiency 




Figure 3: Left : Distribution of NpE as a function of 
^ p. Green and black dots correspond to gamma- 
induced and hadron-induced air shower events, re- 
spectively. Each closed circle with an error bar 
represents the 20%, 50%, 80% of the median dis- 
tribution in each ^ p bin. The solid line shows the 
optimized cut to suppress hadron-induced events. 
Air showers accompanied by no PEs are plotted 
at iVpE = 1.2. Upper right : A'pe distribution in 
the 10 TeV energy band (100 < Y^p < 215). 
Lower right : that in the 100 TeV energy band 
(1000 < Ep < 2154). 

of the 20 inch PMTs used in the simulation are 
shown in Figure 2. After simulating Cherenkov ra- 
diation, propagation of Cherenkov photons in wa- 
ter and the response of the PMTs, the number of 
the collected photoelectrons (A'pe) was obtained 
for each muon detector A'pe resolution for one 
vertical penetrating muon is estimated to be 34 



PEs 



+ 180% 

-18% 



from the simulation. 



Results and discussions 

Figure 3 shows the distribution of ^ iVpE as a 
function of ^ p. J2 ^pe denotes the sum of A'pe 
for muon detectors fired with A'pe > 10. 'Y p — 
1000 corresponds to approximately 100 TeV pri- 
mary gamma-ray energy. ^ A'pe allows for the in- 
fluence of accidental muons, which impinges each 
muon detector at the rate of nearly 1 .6 kHz. In se- 
lecting muon-poor events, we set the optimized cut 
condition as shown in Figure 3. 

Figure 4 shows the survival efficiency of the ak 
shower events after the cut, obtained from the 
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Figure 4: The survival efficiency after the cut. 
Green and black circles represent gamma-induced 
and hadron-induced events, respectively. 

simulation. Around p = 1000, the number 
of hadron-induced events are suppressed down to 
0.01 % or less, while gamma-induced events are re- 
tained by more than 83%. Finally, Figure 5 shows 
the attainable integral flux sensitivity of the Tibet 
ASh-MD array to a point-like gamma-ray source. 
Its 5(7 sensitivity in one calendar year will reach 
7% and -20% Crab above 20 and 100 TeV re- 
spectively, and surpass the existing lACTs above 
20 TeV. Furthermore, it may surpass the next gen- 
eration lACTs above 40 TeV. A further discussion 
is conducted in [7]. The Tibet ASh-MD array will 
contribute to a deeper understanding of the origin 
and the acceleration mechanism of cosmic rays in 
cooperation with other experiments. 
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Figure 5: The attainable integral flux sensitivity of 
the Tibet ASh-MD array to a point-like gamma-ray 
source, together with the sensitivity of HESS and 
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